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Abstract 

An angular analysis of A® — > J/ipA decays is performed using a data sample corre- 
sponding to l.Ofb -1 collected in pp collisions at y/s = 7TeV with the LHCb detector 
at the LHC. A parity violating asymmetry characterising the A® — > J/tpA decay of 
-0.04 ± 0.17 ± 0.07 and a A° b transverse production polarisation of 0.05 ± 0.07 ± 0.02 
are measured, where the first uncertainty is statistical and the second systematic. 
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1 Introduction 



For A® baryons originating from energetic 6-quarks, heavy-quark effective theory (HQET) 
predicts a large fraction of the transverse 6-quark polarisation to be retained after hadroni- 
sation [IJ|2], while the longitudinal polarisation should vanish due to parity conservation in 
strong interactions. For A® baryons produced in e _ e + — > Z° — > bb transitions, a substan- 
tial polarisation is measured [3 45], in agreement with the Z°bb coupling of the Standard 
Model (SM). There is no previous polarisation measurement for A® baryons produced 
at hadron colliders. The transverse polarisation is estimated to be (9(10%) in Ref. [6] 
while Ref. [7] mentions it could be as large as 20%. However, for A baryons produced 
in fixed-target experiments (8 10 , the polarisation was observed to depend strongly on 
the Feynman variable xp = 2 p^/ y/s, being the A longitudinal momentum and y/s the 
collision center-of-mass energy, and to vanish at xf ~ 0. Extrapolating these results and 
taking into account the very small xf ~ 0.02 value for A® produced at the Large Hadron 
Collider (LHC) at y/s = 7TeV, this could imply a polarisation much smaller than 10%. 

In this Letter, we perform an angular analysis of A®—> J/ip (fi + fi~)A(v7T~) decays using 
l.Ofb -1 of pp collision data collected in 2011 with the LHCb detector [11] at the LHC at 
y/s = 7TeV. Owing to the well-measured A—t pn~ decay asymmetry parameter (a a) [12] 
and the known behaviour of the decay of a vector particle into two leptons, the final 
state angular distribution contains sufficient information to measure the A® production 
polarisation and the decay amplitudes 13 . The asymmetry of the A decay (acj) is much 



less precisely measured 12 , however by neglecting possible CP violation effects, which are 
predicted to be very small in the SM 14 15 , a a and — aj can be assumed to be equal. 



Similarly, CP violation effects in A® decays are neglected, and the decay amplitudes of 
the A® and A b are therefore assumed to be equal. Inclusion of charge-conjugated modes 
is henceforth implied. The asymmetry parameter atj, in A® — > J/ipA decays, defined in 
Sec. |2j is calculated in many publications as summarised in Table [T] Most predictions lie 
in the range from —21% to —10% while Ref. 17] obtains a large positive value using HQET. 
Note that the theoretical predictions depend on the calculations of the form-factors and 
experimental input that were available at the time they were made. 

It should be noted that A9 baryons can also be produced in the decay of heavier 



Table 1: Theoretical predictions for the A®— > J/ipA decay asymmetry parameter a^. 
Method Value Reference 



Factorisation 
Factorisation 

Covariant oscillator quark model 
Perturbative QCD 
Factorisation (HQET) 
Light front quark model 



-0.1 
-0.18 
-0.208 
-0.17 to -0.14 
0.777 
-0.204 



16 
17 
18 
19 

s 

20 



1 



b baryons [21], where the polarisation is partially diluted [6]. These strong decays are 
experimentally difficult to distinguish from A® that hadronise directly from a pp collision 
and therefore contribute to the measurement presented in this study. 

A sufficiently large A® polarisation would allow the photon helicity in A® — > A'j and 
A®— > yl*7 decays to be probed (6,22 23 . This observable is sensitive to contributions from 
beyond the SM. 



2 Angular formalism 



The A® spin has not yet been measured but the quark model prediction is spin |. The 
A® — > J/ipA mode is therefore the decay of a spin | particle into a spin 1 and a spin 
| particle. In the helicity formalism, the decay can then be described by four M.\ 1 \ 2 
helicity amplitudes (-M + i , -M_i , M._i_ x and A4 + i >+1 ) where Ai (A2) is the helicity 
of the A (J/ip) particle. The angular distribution of the decay (Jy-) is calculated in 
Ref. [13] and reported in Ref. [24] . It depends on the five angles shown in Fig. Q The 
first angle, 9, is the polar angle of the A momentum in the A® rest-frame with respect 
to n = (p A o x pbeam) /\pa° x Pbeam|, a vector perpendicular to the production plane. The 
second and third angles are 9\ and 0i, the polar and azimuthal angles of the proton in 
the A rest- frame and calculated in the coordinate system defined by Z\ = Pa/\pa\ and 
yi = (n x pa) / \n x Pa\- The remaining angles are 82 and 02, the polar and azimuthal angles 
of the positively-charged muon in the J/ip rest-frame and calculated in the coordinate 
system defined by z 2 = pj/^ /\pj/$ \ and y 2 = (n x pj/^ )/\n X pj/^ \ . The angular distribution 
also depends on the four -MaiA 2 amplitudes, on the parameter, and on the transverse 
polarisation parameter Pb, the projection of the A® polarisation vector on n. 

Assuming that the detector acceptance over (pi and (f> 2 is uniformly distributed, the 






ft' 
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Figure 1: Definition of the five angles used to describe the A®—> J/ijj((i + fi )A(pir ) decay. 
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Table 2: Functions used to describe the angular distributions in three dimensions. 
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/i(a b ,r 0) ri) 




u 1 a 


COS 01, 


COS 9 2 ) 





1 


1 


1 
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a b 
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cos 9 






2 


2r 1 - a b 




COS 6 1 ! 






3 


2r -l 




cos 9 cos 






4 


|(1 -3r ) 


1 


|(3 cos 2 


02 -r 
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|(a 6 - 3r x ) 


n 


^(3 cos 2 


9 2 -l] 


COS # 


6 


-|(o; 6 + ri) 




±(3 cos 2 


9 2 -l] 


cos 9i 


7 


-|(l + r-o) 




|(3 cos 2 


9 2 -r 


cos # cos 9i 



analysis can be simplified by integrating over the two azimuthal angles 




= E M\ M + h,o\ 2 > \M_ h0 \\ \M_^ X \\ l^ + i, +1 | 2 ) 

gi(P b ,aA) hi(cos9,cos9i,cos9 2 ). (1) 
The functions describing the decay only depend on the magnitudes of the 

M 

AiA 2 amplitudes, on P b and a^, and on cos (9, cos#i, and cos 9 2 . Us- 
ing the normalization condition |A / f + i | 2 + |7W_i | 2 + \M_i_ 1 \ 2 + |A^ + i +1 | 2 = 1, 
the ft functions can be written in terms of* the following three parame- 
ters: a b = |A^ + i j0 | 2 - |A^_i >0 | 2 + \M_i_ 1 \ 2 - |A^ + i i+1 | 2 , r = |A^ + i j0 | 2 + |-M_i j0 | 2 and 
r\ = |A^ + i | 2 — |.M_! | 2 . The functions used to describe the angular distributions are 
shown in Table |2j Four parameters (P b , a b , r$ and r\) have to be measured simultane- 
ously from the angular distribution. The a b parameter is the parity violating asymmetry 
characterising the A% — > J/ip A decay. 



3 Detector, trigger and simulation 



The LHCb detector 11 is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < r\ < 5, designed for the study of particles containing b or c quarks. The detector 
includes a high precision tracking system consisting of a silicon-strip vertex detector 
(VELO) surrounding the pp interaction region, a large-area silicon-strip detector located 
upstream of a dipole magnet with a bending power of about 4 Tm, and three stations of 
silicon-strip detectors and straw drift tubes placed downstream. The combined tracking 
system has momentum resolution Ap/p that varies from 0.4% at 5GeV/c to 0.6% at 



3 



100 GeV/c, and impact parameter (IP) resolution of 20 \xm for tracks with high transverse 
momentum. Charged hadrons are identified using two ring-imaging Cherenkov detectors 
(RICH). Photon, electron and hadron candidates are identified by a calorimeter system 
consisting of scintillating-pad and preshower detectors, an electromagnetic calorimeter and 
a hadronic calorimeter. Muons are identified by a system composed of alternating layers 
of iron and multiwire proportional chambers. The trigger 25 consists of a hardware stage, 
based on information from the calorimeter and muon systems, followed by a software stage 
which applies a full event reconstruction. 

The hardware trigger selects events containing a muon with a transverse momentum, 
Pt > 1.48GeV/c or two muons with a product of their pi larger than (1.3GeV/c) 2 . In the 
subsequent software trigger, we require two oppositely-charged muons having an invariant 
mass larger than 2800 MeV/c 2 and originating from the same vertex, or a single muon with 
Pt > 1.3GeV/c being significantly displaced with respect to all the primary pp interaction 
vertices (PVs) in the event, or a single muon with p > lOGeV/c and pt > 1.7GeV/c. 
Finally, we require two oppositely charged muons with an invariant mass within 120 MeV/c 2 
of the nominal J/0 mass [12] forming a common vertex which is significantly displaced 
from the PVs. In the A® — > J/ipA selection described below, we use the muon pairs selected 
by the trigger. 

In the simulation, pp collisions are generated using Pythia 6.4 26 with a specific 



LHCb configuration 127] . Decays of hadronic particles are described by EvtGen 128] 
in which final state radiation is generated using Photos [29]. The interaction of the 
generated particles with the detector and its response are implemented using the Geant4 



toolkit 30 as described in Ref. 31 



4 Signal selection and background rejection 

A first set of loose requirements is applied to select A®— > J/ipA decays. Charged tracks 
are identified as either protons or pions using information provided by the RICH system. 
Candidate A baryons are reconstructed from oppositely-charged proton and pion candidates. 
They are reconstructed either when the A decays within the VELO ("long A"), or when 
the decay occurs outside the VELO acceptance ("downstream A"). The latter category 
increases the acceptance significantly for long-lived A decays. In both cases, the two tracks 
are required to have p > 2 GeV/ c, to be well separated from the PVs and to originate from 
a common vertex. In addition, protons are required to have > 0.5 GeV/c and pions to 
have pt > 0.1 GeV/c. Finally, the invariant mass of the A candidates is required to be 



within 15MeV/c of the nominal A mass 12 . To form J/ip candidates, two oppositely 



charged muons with pr(n) > 0.5 GeV/c are combined and their invariant mass is required 
to be within 80MeV/c 2 of the nominal J/ip mass. To improve the A® mass resolution, the 
muons from the J/ip decay are constrained to come from a common point and to have 
an invariant mass equal to the J/ip mass. Subsequently, A® candidates are formed by 
combining the A and J/ip candidates. We require the A and J/ip candidates to originate 
from a common vertex and to have an invariant mass between 5120 and 6120 MeV/c 2 . 
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Moreover, A® candidates must have their momenta pointing to the associated PV by 
requiring cos#d > 0.99 where 9& is the angle between the A® momentum vector and the 
direction from the PV to the A® vertex. The associated PV is the PV having the smallest 
IP x 2 value, where the IP x 2 is the x 2 difference when the PV is fitted with or without 
the A® four-momentum. 

To reduce the combinatorial background, a multivariate selection based on a boosted 
decision tree (BDT) [32 
candidate: cos#d, the IP x 2 the decay time, the vertex x 2 an d the vertex separation (VS) 
X 2 between the PV and the vertex. Here, the VS x 2 is the difference in x 2 between the 
nominal vertex fit and a vertex fit where the A® is assumed to have zero lifetime. Two 
variables are related to the J/0 candidate: the vertex x 2 an d the invariant mass of the 
two muons. The last variable used in the BDT is the invariant mass of the A candidate. 
The BDT is using simulation for signal and sideband data (M(J/ipA) > 5800 MeV/c 2 ) 
for background in its training. The optimal BDT requirement is found separately for 
downstream and long candidates by maximising the signal significance A^ s ; g / ^/N S i g + iVbkg, 
where iV S jg and iVbkg are the expected signal and background yields in a tight signal region 
around the A® mass. These two yields are estimated using the signal and background 
yields measured in data after the first set of loose requirements and using the BDT 
efficiency measured with the training samples. The BDT selection keeps about 90% of 
the signal while removing about 80% (90%) of the background events for the downstream 
(long) candidates. Less background is rejected in the downstream case due to larger 
contamination from misreconstructed B° — > J/ip K® background decays. Candidates with 
5550 < M(J/ipA) < 5700MeV/c 2 are used for the final analysis. The B° J/ipK° 
background peaks well below the lower end of this mass range. 



33 with eight variables is used. Five variables are related to the A® 



5 Fitting procedure 



An unbinned extended maximum likelihood fit to the mass distribution of the A® candidates 
is performed. The likelihood function is defined as 

N 

' " 1 (2) 



Anass = C -^ r X J] [Y^N^M^A))^ , 



where i runs over the events, j runs over the different signal and background probability 
density functions (PDF), Nj are the yields and Pj the PDFs. The sum of two Crystal Ball 



functions 34 with opposite side tails and common mean and width parameters is used to 
describe the signal mass distribution. The mean and width parameters are left free in the 
fit while the other parameters are taken from the simulated signal sample. The background 
is modelled with a first-order polynomial function. The candidates reconstructed from 
downstream and long A combinations are fitted separately taking into account that the 
resolution is worse for the downstream signal candidates. The results of the fits to the 
mass distributions are shown in Fig. |2j We obtain 5346 ± 96 (5189 ± 95) downstream 



and 1861 ± 49 (761 ± 36) long signal (background) candidates. Using the results of this 
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fit, s Weights (w mass ) are computed by means of the sPlot technique 35 , in order to 
statistically subtract the background in the angular distribution. 

To ensure accurate modelling of the signal, corrections to the j?t and rapidity (y) spectra 
are obtained by comparing the simulation with data by means of the sPlot technique. 
For the A® and A particles, the simulated data is corrected using two-dimensional (pt,v) 
distributions in order to better reproduce the data. These distributions do not depend 
on the polarisation and the decay amplitudes but have an impact on the reconstruction 
acceptance. The same procedure is used on the pion of B° — > J/ipK® decays and is 
subsequently used to calibrate the (pp, y) spectrum of the pion of the A®— > J/ipA decay. 

Since the detector acceptance depends on the three decay angles, the acceptance is 
modelled with a sum of products of Legendre polynomials (Lj) 

/ace ^ , CijkLi(cos8)Lj(cos6i)L k (cos62), (3) 

where i and k are chosen to be even or equal to one. Unbinned maximum likelihood fits 
to the simulated signal candidates are performed, separately for downstream and long 
candidates. The maximum orders of the Legendre polynomials are chosen by comparing 
the fit probability. The requirements i < 5, j < 4, k < 5 and i + j + k < 9 are chosen. 
The results of the fit to the acceptance distributions are shown in Fig. [3] 

We then perform an unbinned likelihood fit to the (cos#, cos#i, cos 62) distribution. 
Each candidate is weighted with u> to t = w mass x u> acc where w mSuSS subtracts the background 



and w acc = 1 / facefcos 6 , cos 9i, cos 62) corrects for the angular acceptance 36 . The sum 
of the tfmass weights over all the events is by construction equal to the signal yield, and 
w tot is normalised in the same way. Since the weighting procedure performs background 
subtraction and corrects for acceptance effects, only the signal PDF has to be included in 




5550 5600 5650 5700 



M(J/yA) [MeV/c 2 ] 




5550 5600 5650 5700 



M(JA(/A) [MeV/c 2 ] 



Figure 2: Mass distribution for the A®—> J/tpA mode for the (left) downstream and (right) long 
candidates. The fitted signal component is shown as a solid blue curve while the background 
component is shown as a dashed red line. 
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Figure 3: Projections of the acceptance function together with the simulated signal data for 
(top) downstream and (bottom) long candidates. 



the fit of the angular distribution. The likelihood is therefore 




where % runs over all events. A simultaneous fit to the angular distributions of the 
downstream and long samples is performed. The ola parameter is fixed to its measured 
value, 0.642 ±0.013 (12). 

The accurate modelling of the acceptance is checked with a similar decay, B°—t J/ipK®. 
Here, the angular distribution is known, and B° mesons are unpolarised. These decays are 
selected in the same way as signal, and the fitting procedure described above is performed. 
Agreement with the expected (cos#, cos#i, cos distribution is obtained. 

6 Results 

The results of the fits to the angular distributions of the weighted A b — >■ J/ip A data are 
shown in Fig. [4j We obtain the following results: P b = 0.06 ± 0.06, a b = 0.00 ± 0.10, 
r = 0.58 ± 0.02 and r\ = —0.58 ± 0.06, where the uncertainties are statistical only. 

The polarisation could be different between A b and A b due to their respective production 
mechanisms. The data are separated according to the A® flavour and fitted using the 
same amplitude parameters but different parameters for the A® and A b polarisations. As 
compatible results are obtained within statistical uncertainties, the polarisations of A® and 
A b baryons are assumed to be equal. 
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Figure 4: Projections of the angular distribution of the background subtracted and acceptance 
corrected A®—> J/tpA data for the (top) downstream and (bottom) long candidates. The fit is 
shown as solid lines. 



A possible bias is investigated by fitting simulated samples generated with sizes 
and parameters close to those measured in data. We generate many samples varying a& 
between —0.25 to 0.25 while keeping r equal to —r\, thus keeping |-M + i +1 | 2 and |.M,i | 2 
equal to zero. We find that the fitting procedure biases all parameters toward negative 
values, slightly for Pf, and tq (~10% of their respective statistical uncertainties) and more 
significantly for a b and T\ (~40% of their respective statistical uncertainties). For Pb and 
7"o, the biases do not change significantly when changing the value of used to generate 
the simulated samples. On the other hand, the biases on and r\ do change, and the 
observed discrepancies are treated as systematic uncertainties. Moreover, the statistical 
uncertainties on the four fit parameters are underestimated: again slightly for Pb and tq 
and significantly, by a factor of ~1.7, for a& and r\. 

We correct the measured values and statistical uncertainties of the four fit parameters. 
The corrected statistical uncertainties are obtained by multiplying the covariance matrix 
with a correction matrix obtained from the study of the simulated samples. This correction 
matrix contains on its diagonal the squares of the widths of the pull distributions of the 
four fit parameters. The remaining entries of this matrix are set to zero as the correlation 
matrix computed with the results of the fits of the generated samples is found to be very 
close to the correlation matrix calculated when fitting the data. 

Finally, the corrected result is P b = 0.05 ± 0.07, a b = -0.04 ± O.f 7, r = 0.57 ± 0.02, 
7*1 = —0.59 ± O.fO, where the uncertainties are statistical only. The corrected statistical 



S 



correlation matrix between the four fit parameters (P&, «{,, r , ri) is 



/l 



V 



0.10 

1 



-0.07 
-0.63 
1 



0.13\ 
0.95 
-0.56 

1 / 



Large correlations are not seen between the polarisation and the amplitude parameters. 
On the other hand, the amplitude parameters are strongly correlated with respect to each 
other, and r x being almost fully correlated. 



7 Systematic uncertainties and significance 

The systematic uncertainty on each measured physics parameter is evaluated by repeating 
the fit to the data varying its input parameters assuming Gaussian distributions and taking 
into account correlations when possible. The systematic uncertainties are summarized in 
Table |3j They are dominated by the uncertainty arising from the acceptance function, the 
calibration of the simulated signal sample and the fit bias. The uncertainty related to 
the acceptance function is obtained by varying the coefficients of the Legendre function 
within their uncertainties and taking into account their correlations. For the calibration of 
our simulated data, the uncertainty is obtained when changing the (pt,v) calibrations of 
the A®, A and pion particles within their uncertainties and obtaining a new acceptance 
function. The function that is used to fit the data does not include the effect of the angular 
resolution. The angular resolution, obtained with simulated samples, is negligible for 9 and 
9 2 . However, it is large, up to ~70%, for small values of Q\. The systematic uncertainty is 
obtained by fitting simulated samples in which the resolution effect is introduced. Effects 
of the deviation from an uniform acceptance in <p\ and 02 assumed in Eq. ([I]) are found to 
be negligible. The simplification to use only one component to describe the background is 
found not to bias the result. Other systematic uncertainties are small or negligible. These 
are related to the signal mass PDF parameters, the background subtraction and cxa- The 
uncertainty related to the background subtraction are obtained when varying the obtained 
result of the mass fit and computing the io maS s weights again. The ola parameter is varied 



within its measurement uncertainties 12 



To compare our results with a prediction on a parameter p, we compute the significance 
with respect to a p test value using a profile along p of the likelihood function, i.e. the 
likelihood value obtained when varying p and minimizing with respect to the other 
parameters. A Monte Carlo integration is performed to include the systematic uncertainties 
in the likelihood profiles. We perform the fit to the data when varying all systematic 
uncertainties and obtain a likelihood profile for each fit of the data. The likelihood profile 
which includes all systematic uncertainties is then the average of all the obtained profiles. 
The significance is defined as S(p = Ptest) — v^0°g £(Ptcst) — l°g£(j»o)), where C(po) is 
the likelihood value of the nominal fit. Significances are given in the concluding Section of 
this Letter. 
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Table 3: Absolute systematic uncertainties on the measured parameters. 



Source 


Pb 


a b 


?"0 


n 


Acceptance 


0.02 


0.04 


0.006 


0.03 


Simulated data calibration 


0.01 


0.04 


0.006 


0.03 


Fit bias 


0.004 


0.04 


0.001 


0.02 


Angular resolution 


0.002 


0.01 


<0.001 


0.005 


Background subtraction 


0.001 


0.006 


0.001 


0.005 


ola 


0.002 


<0.001 


<0.001 


0.01 


Total (quadratic sum) 


0.02 


0.07 


0.01 


0.05 



8 Conclusion 

We have performed an angular analysis of about 7200 A®—> J/if) fi~)A(pii~) decays. 
The A®— >■ J/ipA decay amplitudes are measured for the first time, and the A® production 
polarisation for the first time at a hadron collider. The results are 

P b = 0.05 ±0.07 ±0.02, 
a b = -0.04 ±0.17 ±0.07, 
r = 0.57 ±0.02 ±0.01, 
Tx = -0.59 ±0.10 ±0.05, 

which correspond to the four helicity amplitudes 

\M + i fl \ 2 = -0.01 ± 0.04 ± 0.03, 
\M_i fi \ 2 = 0.58 ±0.06 ±0.03, 
LM_i_i| 2 = 0.49 ±0.05 ±0.02, 

2 ' 

|M + i j+1 | 2 = -0.06 ± 0.04 ± 0.03, 



where the first uncertainty is statistical and the second systematic. The reported polari- 
sation is obtained for the combination of A® and A b decays. More data are required to 
probe any possible difference. 

Our result cannot exclude a transverse polarisation at the order of 10% (6]. However, a 
value of 20% as mentioned in Ref. [7] is disfavoured at the level of 2.7 standard deviations. 

For the A® asymmetry parameter, our result is compatible with the predictions ranging 



from -21% to -10% 16 20 but does not agree with the HQET prediction of 77.7% U] 



at 6.1 standard deviations. 
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